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Here we study the effect of La doping in EuO thin films using SQUID magnetometry, muon spin 
rotation (/xSR), polarized neutron reflectivity (PNR), and density functional theory (DFT). The 
//SR data shows that the Lao. 15 Euo. 85 O is homogeneously magnetically ordered up to its elevated 
Tc- It is concluded that bound magnetic polaron behavior does not explain the increase in Tq 
and an RKKY-like interaction is consistent with the //SR data. The estimation of the magnetic 
moment by DET simulations concurs with the results obtained by PNR, showing a reduction of the 
magnetic moment per Laa^Eui-xO for increasing lanthanum doping. This reduction of the magnetic 
moment is explained by the reduction of the number of Eu-4/ electrons present in all the magnetic 
interactions in EuO films. Einally, we show that an upwards shift of the Eermi energy with La or 
Gd doping gives rise to half-metallicity for doping levels as high as 3.2 %. 

PACS numbers: 75.50.Pp, 75.70.Ak, 75.30.Et, 76.75.+i 


I. INTRODUCTION 

Ferromagnetic semiconductors are attracting great in¬ 
terest for spintronic applicationJIlf^. However, their tech¬ 
nological impact is limited by their low Curie tempera¬ 
ture (Tc). Therefore, it is paramount to find strategies to 
increase their Tq while maintaining their semiconductor¬ 
like character. In the group of ferromagnetic semiconduc¬ 
tors, EuO with a Tq of 69 K has long been considered a 
archetypal model system since the exchange interaction 
is well approximated by the Heisenberg model®^. EuO 
has a magnetic moment of ^ 7 //b por Eu^+ atonP and 
a band gap of 1.12 eV at room temperature (RT)P. One 
of the most remarkable properties of EuO is the Zeeman 
splitting of the Eu-Sd states below 69 K due to the in¬ 
fluence of the Eu-4/ states, causing a conduction band 
splitting of about 0.6 e\®. This high spin splitting energy 
when coupled with low n-dop ing w as proven to transform 
EuO into a half-metal systenP^. Furthermore, the high 
localization of the Eu-4/ electrons, which mediate the 
magnetic interactions, renders it convenient to study the 
exchange interaction phenomena. 

There are several ways of increasing the Tq of EuO. 
The first is by altering the stoichiometry of EuO, e.g. in¬ 
tro duciMoxy gen deficiency into EuO as we have recently 
showrP^l^. The second is by applying strairP^, conse¬ 
quently changing the lattice parameter of the EuO sys¬ 
tem. However, first principles calculations indicate that 
the band gap closes at 5 % of isotropic stress and 6 % bi¬ 
axial straiiP^. The third and most promising mechanism 
to increase the Tq is by doping EuO with other elements. 
In the group of the lanthanide dopants gadolinium has 


shown to increase th e Tc to as much as 170 K for a 4 % 
doping concentratiorP^. 

It was claimed recently by Miyazaki et that a 
lanthanum doping concentration of 10 % increased Tq 
to 200 K, while previous work published by Schmehl et 
reported a maximum increase of the Tq to 118 K for 
1 % lanthanum doping. The latter authors also stated 
that no further increase in Tq occurred for higher doping 
levels. Previous measurements have also showed a Tq 
of 104 K for a 2 % La doping sampl^^. This elevation 
of Tc is believed to be caused by the donation of one 
5d electron by the La atom into the conduction band 
strengthening thereby the magnetic interaction between 
the Eu-Eu next-nearest neighbor atoms 

In a previous publication we studied the effect of n- 
type doping via oxygen deficiency on EuO thin films us¬ 
ing //SRP. We have shown that the results for highly 
oxygen deficient films were incompatible with the bound 
magnetic polarons (BMP)PS1 theory and that the RKKY- 
like interaction is the only known model compatible with 
our datsPIl. 

Here we study the effect of La doping in EuO thin films, 
in particular how the doping affects its magnetic moment 
and how the magnetic order persists above 69 K. This was 
done by probing the underlying physics by //SR, PNR 
and SQUID magnetometry measurements. We have also 
used DFT to draw a theoretical picture of how the La 
doping changes the band structure. We also argue the¬ 
oretically that the increase in doping changes the half- 
metallic behavior of the films. 

The paper is organized as follows: we provide a descrip¬ 
tion of the sample deposition and magnetometry (Sec- 
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tion then use a loca l probe to address the exchange 
mechanism (Section [in[ ). In the second half of the paper 
we use density functional theory (Section |IV[ ) and sup¬ 
porting experimental data (Section |V|) to describe the 
half-metallicity in Laa^Eui_a^O and finish with some con¬ 
clusions on the effect of doping (Section [vT[). 

II. SAMPLE DEPOSITION AND 
MAGNETOMETRY 

The polycrystalline samples were deposited at 400 °C 
using a CEVP magnetron sputtering system with a base 
pressure of 5x10“^ Torr. Co-deposition was performed 
using three targets: a 99.99 % pure EU 2 O 3 , a 99.9 % 
pure Eu and a 99.9 % pure La target. The deposition 
was performed by fixing the EU 2 O 3 fiux and by keep¬ 
ing the total flux from the Eu and La constant. The 
La doping was controlled by changing the relative flux 
between Eu and La. The growth was performed in an 
Ar+ plasma at a pressure of 2 mTorr with a flow rate of 
14 seem. The substrates used were Si(OOl). A Si 02 layer 
was deposited between the substrate and the La^^Eui-a^O 
film and another on the top as a capping layer to prevent 
further oxidation of the film. We deposited 10 %, 15 % 
and 20 % La doped EuO samples. 

The SQUID magnetometry measurements were per¬ 
formed in a Quantum Design MPMS 2 system. Each 
sample was cooled in a field of 50 Oe to 5 K followed 
by the M{H) measurements at increasing temperatures. 
Einally each sample was cooled at zero field and the tem¬ 
perature dependent measurements where performed at 
different fields. 

Eigurejl] shows the SQUID data for the 15 % La doped 
EuO sample. The hysteresis curve shows a coercive field 
of 75 Oe compared to 54 Oe for the pristine EuO sample 
reported in Ref. [H The inset shows the bulk moment ob¬ 
tained at different saturation fields. The inset of Eigure 
also shows the development of a clear ferromagnetic 
phase below 96 K, with a paramagnetic response above 
this temperature.The decrease of the coercivity with tem¬ 
perature is shown in Eigure 

Similar measurements of the 10 % and 20 % samples 
find a Tq of about 94 K and 84 K, respectively. The lower 
value for the 20 % La doped sample is due to disorder 
and the formation of impurities. The paramagnetic sus¬ 
ceptibility is well described by the Curie-Weiss law and 
we used this to determine Tq. 


III. MAGNETIC SPACIAL HOMOGENEITY 

The elevated Curie temperature of the EuO system 
has attracted considerable attention with s evera l mod¬ 
els being invoked to describe this behavio j^^ l ^^ i Low 
energy muon measurements have been shown to be a 
highly effective local probe of the magnetism of such sys¬ 
tem^. Our measurements were performed at the Paul 
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FIG. 1. Temperature dependent SQUID data for the 15 % La 
doped EuO sample. The inset shows the temperature depen¬ 
dence of the magnetic moment. The moment was normalized 
by the PNR results at 5 K. 



FIG. 2. Goercivity as a function of temperature for the 15 % 
La doped EuO sample. 


Scherrer Institute, Switzerland, in zero applied field (ZE) 
and weak transverse field (wTE = 28.2 G) configurations, 
where muons were implanted into the film at various en¬ 
ergies between 6 and 14 keV. In /iSR experiments^ spin 
polarized positive muons are implanted into the sample, 
where they stop rapidly at interstitial sites of high elec¬ 
tron density, and their spin direction evolves in the mag¬ 
netic field at their stopping site. Each implanted muon 
decays with a lifetime of 2.2 /is emitting a positron pref- 
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FIG. 3. Raw muon data for the Lao. 15 Euo. 85 O sample: (a) 
Zero-field measurements with the fits to Eq. (b) Weak 
transverse field measurements with the fits to Eq. 


erentially in the direction of its spin at the time of de¬ 
cay. The evolution of their spin polarization as a func¬ 
tion of time is determined by measuring the direction 
of the emitted positrons. In low-energy muon experi- 
ment^^^H^, accelerating electric fields are used to con¬ 
trol the implantation depth of the muons after they have 
passed through a cryogenic moderator. Example data are 
shown in Fig. The lack of clear oscillations below Tq 
in ZF, compared with the data for bulk EuO obtained by 
Blundell et may be due to small variations in the 

internal fields caused by a distribution of grain bound¬ 
aries and defects in the film. This was also evident in our 
previous measurements of EuOi_a;^. 

We described the ZF data using the functioiP: 

A{t) = Ai exp(-At) + A 2 exp(-At) + ^bg, (1) 

where the sum of Ai and A 2 (17.5 %), and the back¬ 
ground contribution A\yg (6 %) are fixed by the geometry 
of the spectrometer and sample. The first term, describ¬ 
ing the fast relaxation A, is found to be zero above Tq 
and represents an incoherent precession of muons about 
fields perpendicular to their spin polarization. This re¬ 
laxation rate A, shown in Figure (a), will therefore 
vary with the size of those fields and gives us some in¬ 
sight into the magnetic order parameter, albeit far more 
limited than a well-defined oscillation frequency. This 
change in A exhibits a similar temperature dependence 
to the magnetization recorded in our bulk measurements, 
fully consistent with following a ferromagnetic order pa¬ 
rameter by both techniques. The second term describes 
a slow relaxation of muon spins. A, which below Tq we 
can attribute to fluctuations of magnetic fields parallel 
to the muon spin polarization and above Tq is caused by 
paramagnetic spin fluctuations. Figure |^b) shows that 
there is a small peak in A around Tq in Lao. 15 Euo. 85 O 
that is less pronounced than that in pristine EuOEI. 

Below Tq we note that A 1 /A 2 ^ 2 , which is the value 
anticipated for a polycrystalline ordered magnet. Us¬ 
ing this we can estimate the fraction of the probed sam¬ 
ple volume entering a magnetically ordered state at each 
temperature as Pmag = 1-5 x Ai/{Ai + A 2 ). In Fig¬ 
ure. (c) the static magnetic volume Pmag is shown as 
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FIG. 4. Fitted parameters for the zero field muon data (a- 
c) and weak transverse field data (d-f): (a) Fast relaxation 
rate, A. (b) Slow relaxation rate, A. (c) Static magnetic vol¬ 
ume fraction, Pmag- (d) Magnetic field, B, experienced by the 
implanted muons, (e) Relaxation rate, 77. (Square-root expo¬ 
nential for Lao. 15 Euo. 85 O and simple exponential for EuO.) 
(f) Relaxing asymmetry, Ar. The vertical dotted lines denote 
the Gurie temperatures for EuO (68 K) and Lao. 15 Euo. 85 O 
(96 K), with data for EuO taken from Ref. O 


a function of temperature. These data shows that qua¬ 
sistatic magnetic fields develop through the whole sample 
volume, within the ^ 5 % measurement resolution, at Tq- 
The wTF measurements at low-temperature show a 
fast relaxation and a slowly-relaxing oscillation, which 
can be attributed respectively to muons experiencing the 
large spontaneous fields and the weak applied field. This 
gives us another probe of the volume of the sample in 
which the muons are implanted that is magnetically or¬ 
dered. To simplify the determination of this volume frac¬ 
tion we fitted the data omitting the first 0.25 /is to the 
function: 


A{t) = Ar exp(-v^) cos( 7 ^P + 0). ( 2 ) 

This describes the slow relaxing precession of muons not 
experiencing large internal magnetic fields within the 
magnetic volume of the sample. The parameters 77 and 
B are shown in Fig. (d) and (e), respectively. The be¬ 
havior of all three parameters is broadly similar to the 
equivalent parameters for EuO, except that transverse 
field relaxation here takes a square root exponential form, 
rather than the simple exponential observed for EuOi_a^. 
This difference can be explained by the presence of an in¬ 
homogeneous distribution of fields at muon stopping sites 
due to non-magnetic site dilution by the La doping. The 
amplitude A^ [Fig. |^f)] drops on entering the magneti¬ 
cally ordered phase, with the remaining amplitude being 
consistent with the background determined for the ZF 
measurements. 
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Both ZF and wTF data show that the sample suffers 
a sharp transition from a fully polarized state to a para¬ 
magnetic one. This behavior is similar to what we pre¬ 
viously obtained for EuOi_J^. Therefore, we conclude 
that BMP could not be present in Lao. 15 Euo. 85 O and the 
increase of Tq is a consequence of the existence of the 
RKKY-like interaction between the 4/ and the bd states. 


IV. HALF-METALLICITY IN La^Eui_,,0 

Of central interest in the EuO system is the half- 
metallic behavior below the Curie temperature. To inves¬ 
tigate this as a function of La doping density functional 
theory calculations were performed using the CASTEP 
plane-wave pseudopotential packag^^. As noted in our 
previous worlPI, bulk EuO is poorly described by stan¬ 
dard DET employing local or semi-local functionals, 
which describe it as a paramagnetic metal instead of 
a ferromagnetic insulator. We use instead the DET+U 
formalisnP^, which augments the Local Spin Density Ap¬ 
proximation (LSDA) with a Hubbard model description 
penalizing non-integer occupation of the Eu 4/ orbitals. 
As in the previous work, we used U = 7.3 eV as it accu¬ 
rately reproduces both the experimental bandgap and 
lattice parameter of bulk EuO. Simulations were per¬ 
formed with a plane-wave basis with a cutoff energy of 
600 eV, and a Monkhorst-Pack A;-point mesh with a min¬ 
imum spacing of 0.1 A“^. 

We first constructed a 64-atom supercell comprising 
2 x 2 x 2 copies of the 8 -atom conventional cubic cell. 
Within this cell the calculated magnetic moment per Eu 
atom is exactly 7.0 /ie- La doping was then introduced 
by substituting a chosen number of Eu atoms with La 
atoms and relaxing the resulting geometry. 

Since La does not itself have occupied 4/ states, each 
La substitution might be expected to reduce the overall 
magnetic moment by 7.0 /ie- However, in fact the mag¬ 
nitude of the reduction will be less than this, as the La 
acts as an electron donor to the bd states comprising the 
conduction band. The lowest-lying conduction band (at 
F) is spin-split in the simulations by 0.7 eV due to the 
effective field from the aligned Eu 4/ moments. There¬ 
fore the electrons occupying the conduction band are at 
low concentration fully aligned with the magnetic mo¬ 
ment of the 4/ electrons. As the concentration increases, 
however, the minority-spin channel begins to be filled as 
well, and the contribution per conduction electron falls. 
At this point the half-metallicity is lost as the Eermi level 
crosses both majority and minority spin states Eig.[^(c). 

Eigure shows the effective band structure (EBS) of 
the La-doped super cell calculations, projected into the 
Brillouin zone of the conventional cubic unit cell of per¬ 
fect EuO. We use the method of Popescu and Zungei!^ 
for the EBS plots, as implemented by Brommer et al^. 
As the La-doping fraction x is increased, we see the Eermi 
level rise further into the conduction band around F. Si¬ 
multaneously, the 4/ bands become less distinct, as the 



FIG. 5. Effective band structure for a 64-atom EuO supercell: 
(a) Pristine EuO. (b) 1 La substitutional atom (3.125%). (c) 
2 La substitutional atoms (6.25%). (d) 3 La substitutional 
atoms (9.375%).The La substitutional dopants on Eu sites 
are projected onto the 8-atom EuO conventional cell. The 
left panel corresponds the majority-spin states and the right 
the minority-spin states. The blue line is the Fermi level. 
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FIG. 6. Effective band structure for a 64-atom EuO supercell 
with 1 Gd substitutional atom (3.125%). The substitutional 
dopants on Eu sites are projected onto the 8-atom EuO con¬ 
ventional cell. The left panel corresponds the majority-spin 
states and the right the minority-spin states. The blue line is 
the Eermi level. 


La dopant atoms both disrupt the periodicity of the sys¬ 
tem and reduce its overall magnetic moment. At high 
X, we see evidence of subbands forming in many regions 
of the plot, but this is due to the artificial periodicity of 
the 64-atom cell and would not be observed so clearly for 
truly random distributions of defects in larger cells. 


Mulliken population analysis using pseudo-atomic or¬ 
bitals as projectors indicates that the magnetic moment 
is not strongly localized on La sites. These only have a 
net spin in the range 0 . 1 - 0.2 /ip depending on the number 
of La atoms present in the cell. Instead, the net magnetic 
moment is distributed over the whole cell, indicating that 
it occupies rather delocalized bd orbitals. 


Given the high ordering temperature with Gd doping 
it is also opportune to compare the effects of La and Gd 
doping in order to better understand the phenomena of 
electron doping. We have therefore performed DFT cal¬ 
culations substituting a number of Eu atoms for Gd, and 
we show effective band structures from these simulations 
in Fig. In comparison to the La-doped simulations, 
doping with Gd shows a overall similar behavior in re¬ 
lation to the position of the Fermi level with respect to 
the majority and minority spin bands, but the 4/ density 
of states remains constant for increasing doping levels as 
there is no decrease in 4/ occupancy per unit cell. There¬ 
fore, appropriate levels of Gd doping can also be expected 
to result in a half-metallic material. 



FIG. 7. PNR data and analysis for the pristine EuO, 10 % 
and 15 % La doped EuO samples at 5 K. The spin-up and 
spin-down data is shown in red and blue, respectively. 


V. VARIATION OF THE MAGNETIC MOMENT 
WITH La DOPING 


To investigate the critical effect of the La doping in 
EuO thin films one needs to understand the impact it 
has on the overall magnetic moment and also if any sub¬ 
stantial changes occur in the bandstructure. The ex¬ 
perimental determination of the magnetization of the 
Laa,Eui_a,0 films was done using low temperature PNR 
on the neutron refiectometer CRISP at ISIS, UK, follow¬ 
ing preliminary measurements on beamline D17 at the 
ILL, France. The data analysis was performed by us¬ 
ing the GenXpSI XPolly fitting softwares. Figure 
shows the PNR reflectivity obtained for the three mea¬ 
sured Laa,Eui_a,0 samples at 5 K, showing that the 15 % 
La doping sample has a reduced magnetic moment of 5.7 
/iB per Lao. 15 Euo. 85 O formula unit in comparison with 
the 6.98 /iB per stoichiometric EuO formula unit. 

Figure compares the magnetic moment calculated 
using DFT and measured by PNR. The figure shows a 
good agreement between the theory and the experiment 
at 10 % doping. The 15 % lanthanum doped EuO sam¬ 
ple has a magnetic moment slightly reduced compared 
with the DFT calculation and we believe this is owing 
to the fact that the random distribution of the La atoms 
creates disorder that weakens the magnetic interactions. 
Another explanation is the presence of a small percent¬ 
age of impurities lowering the overall magnetic moment 
measured by PNR. The figure also shows the calculated 
maximum and minimum possible polarization of the con¬ 
duction band, in a model which assumes that at concen¬ 
tration X of La substituting Eu, the moment of the 4/ 
electrons will be given by 7(1 — x) /ib- The maximum 
polarization curve assumes a further x jj^B contribution 
from fully-polarized electrons in the conduction band. 















6 



FIG. 8. Comparison between the DFT calculations and the 
PNR analysis. The blue-dashed and red-dotted-dashed lines 
corresponds to the maximum and minimum polarization com¬ 
ing from the doping electrons, respectively. 

while the minimum polarization curve assumes they are 
fully unpolarized. As explained before, the doping elec¬ 
trons populate only the majority-spin states at low dop¬ 
ing levels. However as the doping increases, the Fermi 
energy is shifted upwards, crossing the minimum of the 
minority-spin conduction bands. Beyond this point both 
minority and majority spin states will be occupied and 
the overall polarization drops. This can be seen in Fig¬ 
ure the DFT and maximum polarization curves start 
at similar values, and at about x = 0.0625 they begin 
to diverge. This divergence is attributed the loss of the 
half-metallicity that occurs ai x = 0.03125. 

VI. CONCLUSIONS 

In this paper we have shown that the La doping acts 
like an n-dopant increasing the density of carriers in the 
bd conduction band and strengthening the Af-bd (Eu-Eu) 
interaction. The /iSR experiments show that the elec¬ 
tron doping in Laa,Eui_a,0 thin films has a similar effect 
on the magnetic properties to that which was found for 
EuOi-a,. One of the indicators is the fact that 
develops an abrupt transition at the elevated Tq (96 K). 
This is a crucial indicator that there is no magnetic phase 


separation by the presence of bound magnetic polarons. 
It also indicates that no significant clustering of La is 
occurring, as this would broaden the transition. In the 
absence of bound magnetic polarons the electrons pop¬ 
ulate the bd band and the RKKY-like interaction is the 
one that explains the Tq enhancement at high doping 
levels. 

The PNR and DET show a reduction of the overall 
magnetic moment of the Laa,Eui_a,0 formula unit. This 
is attributed to the replacement of Eu by La atoms re¬ 
ducing the total amount of 4/ electrons in the system. 
The 4/ electrons are responsible for carrying the mag¬ 
netic moment and when removed a dilution of the total 
magnetic moment in the system occurs. The difference 
between the PNR and DET at higher doping levels is 
ascribed to disorder and impurities in the sample. 

The DET calculations predicts that the band structure 
remains almost unaffected for increasing doping where 
the only change is an upwards shift of the Eermi en¬ 
ergy. Eor doping levels below 3.125 % this change is just 
enough to make the Eermi energy intersect the majority 
spin states leaving the minority spin stated unoccupied 
- forming a fully spin polarized conduction band. This 
half-metallic behavior occurs at much higher levels than 
the ones reported for oxygen deficient EuO, where at the 
same vacancy doping levels the system was already show¬ 
ing both the majority and minority channels occupied. 
The DET calculations also show that the range of dop¬ 
ing where half-metallicity can be expected in La or Gd 
doped EuO is limited by the size of the conduction band 
spin splitting. The presence of half-metallicity with La or 
Gd doping makes EuO a extraordinary system to study 
the phenomena of half-metallicity in strongly correlated 
systems. 
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